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Abstract 

O 

A cosmological model with a cyclic interpretation is introduced, which is subject 
5h \ to quantum back-reaction and yet can be treated rather completely by physical co- 

herent state as well as effective constraint techniques. By this comparison, the role 
of quantum back-reaction in quantum cosmology is unambiguously demonstrated. 
Also the complementary nature of strengths and weaknesses of the two procedures 



wo: 



> 

is illustrated. Finally, effective constraint techniques are applied to a more realistic 
model filled with radiation, where physical coherent states are not available. 

ON 

1 Introduction 

To understand the generic behavior of a quantum system, effective equations^] are useful. In 
contrast to individual wave functions, or even just stationary states, they directly provide 
(approximate) equations for time-dependent expectation values. Since the dynamics of 
expectation values depends on the whole motion of a state — by quantum back-reaction, 
all moments of a state couple to the expectation values - - these equations in general 
differ from the classical ones by quantum corrections. If quantum corrections, for instance 
an effective potential, can be found explicitly, an interpretation of quantum dynamics in 
generic terms becomes much easier. Such results are more general compared to conclusions 
based on individual states. 

Especially in quantum cosmology, the ability to draw generic conclusions is important. 
Not much is known about the state of the universe except, perhaps, that it currently 
can well be considered semiclassical. But semiclassicality is not a sharp notion, and so 
wide classes of states, differing for instance by the sizes of their quantum fluctuations 
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1 We use the term "effective equations" in a general sense to encompass equations that might follow 
from an effective action, or be derived in a canonical way. Systematic procedures exist for both ways, and 
the results agree in cases where they have been explicitly applied [I] [2] • 
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or correlations, are still allowed. In any such situation, a generic analysis is called for, 
most crucially when long-term evolution is involved, or when one evolves toward strong- 
curvature regimes such as the big bang singularity where quantum effects of all kinds are 
expected to be important. 

It is sometimes suggested, at least implicitly (and especially in the context of loop 
quantization), that quantum cosmology might somehow be different from other quan- 
tum systems, and that quantum back-reaction could be ignored in its effective equations. 
Quantum back-reaction might be weak for certain states or in certain regimes, especially 
for models close to solvable ones, but this observation cannot be generalized. Like the 
harmonic oscillator in quantum mechanics, there are harmonic cosmologies [3l H] where 
expectation values of states follow exactly the trajectories of a corresponding classical sys- 
tem. Such systems are entirely free of quantum back-reaction. For "small anharmonicity" , 
quantum back-reaction might still be weak, as it is realized in quantum cosmology for mat- 
ter dominated by its kinetic energy density [5] . But the tough reality of stronger deviations 
from the solvable ideal of harmonic systems can introduce severe quantum back-reaction, 
which must be studied in an unbiased and systematic way. 

Here, we introduce a model of quantum cosmology which is not solvable but still treat- 
able by two rather different methods: effective constraints and physical states. The model is 
an anisotropic cosmology of locally rotationally symmetric (LRS) Bianchi I symmetry type, 
filled with an isotropic, slightly sub-stiff fluid of negative energy density p(a) oc — log(a) 2 /a 6 
where a is the average scale factor of the anisotropic geometry. With this specific density, 
the model becomes treatable by physical coherent states, which justifies its contrived and 
exotic form. Effective constraint techniques [61 [7j, applicable much more widely, do not 
require such a tailored matter source; they are considerably more powerful. As we will 
see by the explicit comparison of this paper, they capture the information in semiclassical 
physical states to an excellent degree. Moreover, effective techniques self-consistently de- 
termine their ranges of validity. The general applicability of effective constraints will be 
demonstrated by our final analysis of a model whose matter content, more realistically, is 
pure radiation. 

2 The model 

An anisotropic Bianchi I model has a line element 

d S 2 = -di 2 + ^a^) 2 (dx J ) 2 

7=1 

with three independent scale factors aj(t) as functions of proper time t. For a canonical 
formulation, we denote the momenta of hj := a 2 by ir J . It is convenient to introduce Misner 
variables (a, /?+,/?_) and their momenta (p a ,p + ,p-) by the canonical transformation [8] 

I III 

-log(7ii/ajj) =: a + (3 + + VS(3_ , 2h X 'K 1 = : -p a + -p + + — =p_ (1) 
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- log(/i 2 / a 2 Q ) =:a + /3+- V3f3_ , 2/i 2 vr 2 =: ^p a + -p + - -^j=P- ( 2 ) 
1 11 

- log(/i 3 /ao) = : « - , 2/i 37r 3 =: -p a - ~p+ ■ (3) 



In these definitions, ao is a reference scale factor (e.g. = yhji^hz at one moment of 
time) introduced to be insensitive to coordinate rescalings. 

Canonical dynamics in general relativity is determined by the Hamiltonian constraint 



^ ab -\^\? y/deU (3) 
Vdet/i (167rG) 2 



with the spatial metric /i a ft, its Ricci scalar and its momenta 7r afc . Reduced to Bianchi I 
metrics in Misner variables, specified to a lapse function N = Vdet h = aia 2 a^, it simplifies 
considerably to the form 

C B ianchiI = 24 (-2^ + P+ + P-) ■ ( 5 ) 

As one can check directly, the constraint generates the correct Hamiltonian equations of 
motion in coordinate time, from which the Kasner solutions follow. 

We now restrict the model further by requiring the anisotropy parameter (3- and its 
momentum p_ to vanish: (3- = = p_. In this way, which can easily be confirmed 
to be consistent with the equations of motion, we enhance the symmetry and leave two 
independent gravitational variables: the logarithm of the average scale factor, a, and one 
anisotropy parameter j3 + . The resulting Hamiltonian constraint is equivalent to that of a 
free, massless relativistic particle. 

To introduce a "potential" , we will work with a matter source whose energy density is 

log(a/a ) 2 6 _ 6q 2 / r n 

(a/a ) b 

where, in addition to a Q already introduced, p is the matter energy density at some 
"initial" time. In the presence of matter, its density multiplied with a 6 is to be added to 
the constraint §5§. It become^] 

C L RS = -Pa + P\ ~ a 2v ( 7 ) 

if we make convenient (and irrelevant) choices for the prefactors and define v = po a o- O ur 
constraint then becomes that of a "relativistic harmonic oscillator" as studied in [7]: 

pl=pl + ua 2 . (8) 



2 If we introduce a finite region of coordinate size Vo to integrate out homogeneous quantities such as 
the symplectic form f d 3 xSh a b A STr ab = Vodhi A Six 1 = Shi A Sw 1 , the momenta tt 1 — Vqtt 1 depend on the 
rather arbitrary Vo- The lapse function as chosen here would then be homogeneous in Vq, too, such that 
V^a G p is the matter contribution to the constraint. With this, all its terms scale in the same way if Vo is 
changed. Solving the constrained system, we obtain the same reduced phase space for all choices of Vq. 
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In this analogy, we consider (3 + as our time variable, and p + as the corresponding "energy." 
Evolution of a and p a with respect to f3 + is then generated by the Hamiltonian p + = 

A more realistic matter content could be chosen as radiation, with an energy density 
p oc a" 4 = e~ Aa . Here, the Hamiltonian is p + = ±^Jp\ — ue 2a . 



2.1 Classical behavior 

We define H = \/p 2 a + va 2 , such that our constraint solved for p + takes the deparameter- 
ized form 

Cdcparamctrizcd = P+ =t H(a,p a ) =0. 

(From now on we choose p + positive, to be specific.) This constraint generates Hamiltonian 
equations of motion for the anisotropies p + = and $ + — 1, allowing us to identify the 
time parameter along its flow with /3 + as an internal time. All derivatives in Hamiltonian 
equations of motion are then with respect to f3 + , specifically 

{a,H} = — , —— = {p a , H} = — 



d/3+ H ' dp+ J if 

Since H = is constant in we can combine these equations to a second order one 
for a([3 + ), d 2 a/d/3l = —ua/H 2 solved by 

a(P+) = A sm((3 +v ^/H) + B cos((3 +v ^/H) (9) 

with integration constants A and B. The equation of motion for a + tells us that 

p a {(3 + ) = {Acos(P + V^/H) - Bsmtf+y/V/H)) . (10) 

With this, the integration constants can be related to H by H 2 = v (A 2 + B 2 ). Solutions 
in phase space are ellipses of axis lengths H = \p + \ and y/vH, traversed in time f3 + with 
frequency {2nH / yju)^ 1 . 

To derive the behavior with respect to proper time, we use the original constraint Clrs 
not yet deparameterized, and remember that we chose a lapse function iV = 010203 = a\e Za . 
Thus, the equation of motion for f3 + with respect to proper time r is 

d P+ _ f o c -3a n \ _ o -3« -3 

-^T-|P+,e G L Rsj-^e p + a . 
From here, we determine proper time as a function of (3+ by integrating 



2p 



+ 



-,3a 



^d[3 + (11) 



with our solution a{(5 + ). Inverting the solution allows us to insert P+(r) into cx(/3+) and 
p a (/3 + ), resulting in solutions as functions of proper time. 
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It is not easy to integrate t(/3+) explicitly, but it is clear from the integrand that 
t(/3 + ) is a monotonic, one-to-one function which can be inverted globally. Thus, ot(r) 
and p q (t) are defined and finite for all values of proper time. There is no singularity 
in this model. Clearly, the negative amount of matter energy violates energy conditions 
sufficiently strongly to provide the cyclic bouncing solutions embodied by ellipses in the 
(a, p a ) -plane. 

One can see this directly from the Friedmann-type equation resulting from the Hamil- 
tonian constraint. We have p 2 a = p 2 + — ua 2 , where p + is constant and p a can be obtained 
from the equation da/dr = {a, 3 e~ 3a CLRs} = — 2<2q 3 e~ 3a p a , relating it to the derivative 
of a by proper time r. Thus, the constraint equation reads 



and for a = a exp(a) implies 



A 2 A pI a (loga/ao) 2 no s 
4 Po 7Z TZ \ fi • \ lz ) 



a J a 6 (a/aoY 

On the right hand side, we identify Ap 2 + /cP as the anisotropic shear term, and the contri- 
bution — 4p ag (log a/a ) 2 /a 6 as our exotic energy density. 

The Friedmann equation shows when a can vanish, which is realized for loga/ao = 
±p + / y/u, two solutions which correspond to the maximal and minimal a along our cir- 
cles. The extrema indeed give us a maximum for loga/ao = p+j\fv and a minimum for 
log a/a = —p+j'sjv: The sign of a, by the Raychaudhuri equation, is determined by the 
sign of p + 3P of all the matter sources combined. Here, pressure is obtained from the 
energy density as the negative derivative of energy by volume, i.e. by 

_dE d(a 3 p(a)) 
dV 3a 2 da 1 J 

resulting in P neg = p neg (l — f (loga/ao)" 1 ) for the negative energy fluid with p neg = 
— 4i/(log a/a ) 2 /a 6 , and P s hear = Pshear for the shear contribution with p s h ea r = 4p^/a 6 . 
For p + 3P, this gives 



P + 3P = p ncg 4 - - — — — r + 4p shear 
V. log(a/a )/ 

which at the extrema log a/ao = ±p+ / \fv evaluates to ±8y/vp+ / a 6 . The extrema provide a 
maximum at log a/a = p + /y/u and a minimum at log a/a = —p + /y/D ) and the evolution 
is that of a cyclic model with infinitely many bounces and recollapses. 

For the model with radiation, we have equations of motion da/d/? + = p a /H and 
dp a /d(3 + = ve 2a /H, with H still constant. The second order equation for a{(3 + ) becomes 
d 2 a/d(3 2 _ = ue 2a /H 2 . For the solution of this model, see Section HI 
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2.2 Quantum representation 

To represent the model as a quantum system, we start with the obvious kinematical Hilbert 
space L 2 (M 2 , dad/3+) of square integrable wave functions of two variables, /3+ and a. The 
momentum operators are derivatives times — iH, as usual. To arrive at observable informa- 
tion and physical states, we have to implement the constraint operatoiH 

C = pl-pl-& 2 = " ft2 g^ + ft2 g^5 " " 2 > ("J 

find its kernel and equip it with a physical inner product. 

To analyze the quantum constraint, we reformulate it in a first-order way in our time 
variable f3 + by taking a square root: 

- ?^-*(a, P+) = ±{f a + « 2 ) 1 *(a, (3+) =: ±tf¥(a, /?+) , (15) 

introducing the deparametrized Hamiltonian operator H = (p 2 a + a 2 ) 1 / 2 . All solutions of 
this equation can be expressed as 

oo 

*±(a,/3+) = ^c n v9 n (a)exp(=FiA n /?+//i) (16) 

n=0 

where y? n are the eigenstates of H with eigenvalues h\ n . Since if is the square root of the 
(positive definite) harmonic oscillator Hamiltonian (with "mass" m = 1/2 and "frequency" 
u — 2), its eigenstates are of the well-known form, with eigenvalues X n = \J (2n + Tjh. 

In f|T6|) . the subscript ± indicates the sign taken when solving for p + by a square root. 
Solutions naturally split into two classes, positive-frequency solutions and negative- 
frequency solutions The separation into two classes becomes relevant when we intro- 
duce the inner product 

(M/,$) aux :=z J da (**( a ,P + )^L$(a,P + ) -$(«,/?+) (17) 

of Klein-Gordon form. Although (3 + appears on the right hand side, the inner product 
evaluated on solutions (fT6"j) is time independent. However, (-, -) aux is not positive definite: 
it is positive for positive-frequency solutions, but negative for negative-frequency solutions. 
A positive-frequency solution is automatically orthogonal to a solution with the sign of its 
frequency flipped. To correct the sign, we define the physical inner product as 

(\]>, $) aux if \1> and <E> are positive-frequency 
(vp, $) := ^ — (vf, $)aux if ^ and $ are negative-frequency (18) 
otherwise 



3 Henceforth, for this model, we drop v for simplicity. It can be absorbed into the variables by first 
dividing the constraint throughout by y/V, followed by the canonical transformation a' = i/ia, p' a = v~*p a , 
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and extend it linearly to superpositions of positive and negative frequency solutions. (Al- 
ternatively, we may declare positive and negative frequency solutions, respectively, to define 
two superselection sectors. The procedure here is analogous to [9].) This completes the 
construction of the physical Hilbert space. 

From physical states, /3+-dependent expectation values (or moments) can be computed, 
playing the role of evolving observables. As with the non-relativistic harmonic oscillator, 
this is most easily done using ladder operators: a = y/h/2 (a) + d) , p a = i^/h/2 (a) — d) . 
A direct calculation gives 

oo 

= ^ 2h ( n + X ) Re (CnCn+l eX P (-«/?+( A n+ i - \n)/h)) (19) 

n=0 

oo 

(# + ,p a * + )(/3+) = V^Kn + l)Im (cnCn+i exp (-i(3+(\ n+1 - X n )/h)) (20) 

n=0 

and similarly for moments. 

Physical Hilbert spaces can be constructed and physical states decomposed in this way 
whenever one knows an explicit diagonalization of the Hamiltonian H. For our model, 
the closeness to the harmonic oscillator has an additional advantage in that it allows 
us to use its simple form of coherent states — Gaussians of arbitrary width expanded 
in the stationary states — as initial values for evolution in /3+. For the non-relativistic 
harmonic oscillator, the resulting physical states would be dynamical coherent states: their 
shape would remain unchanged and they keep saturating the uncertainty relation at all 
times. Moreover, their expectation values follow the classical trajectories exactly, without 
quantum back-reaction. 

For the relativistic harmonic oscillator, and thus our anisotropic toy model, the dy- 
namical behavior of the states is still to be seen. We thus assume an initial state, at some 
fixed value of j3 + , of the kinematical coherent form: 

c n = exp -i-L. _ , z G C (21) 



such that ty(a, 0) = (2/tt) 1/4 exp (-|(M 2 - z 2 + 2a 2 - Aiza)) from (ilSl) . A time-dependent 
physical state then has coefficients 



c e -iKP+/n = _J_ e ~\z\ 2 /2 z n e -iV2E+l/3+/Vh 



in its expansion by the (p n (a). With the square root of 2n + 1, the exponentials cannot 
simply be combined to a /5 + -dependent z(f3 + ) n . The shape of the state changes as (3 + moves 
away from zero: the time-dependent coefficients are no longer of the form Q2ip for (3 + ^ 0. 
Physical states with initial conditions given by the coherent states of the non-relativistic 
harmonic oscillator are not dynamical coherent states. The model introduced here does 
show spreading and quantum back-reaction, which makes it interesting for a comparison 
with effective constraint techniques. 
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2.3 Effective constraints 



In an effective treatment of a quantum system, we can consider the same dynamics, but 
focus on the algebra of observables. Results will thus be manifestly representation indepen- 
dent. To set up this framework, no approximations are required; we are thus dealing with 
an exact quantum theory. Only when evaluating the equations, which would give us expres- 
sions such as (a)(/3 + ) or moments of a state as functions of (3 + , do approximation schemes 
typically enter. This is no difference to a representation dependent treatment, where exact 
evaluations of expectation values such as ( 1191) or ( 1201) are hard to sum explicitly. 

At the kinematical level, effective techniques are based entirely on the algebra of basic 
operators, in our case [/3 + ,p + ] = ih and [a,p a ] = ih, with all other basic commutators 
vanishing. As it happens at the quantum level, dynamics is brought in by a constraint 
operator C which might have more complicated algebraic relationships with the basic 
operators, no longer forming a closed algebra. 

A whole representation of these algebraic relationships on wave functions carries much 
more details than necessary for extracting physical results. Instead, it is often convenient to 
focus directly on expectation values and derive dynamical equations for them, avoiding the 
detour of computing wave functions. Expectation values are not sufficient to characterize 
a state or its dynamics, but when combined with all moments 

A(Oi...O n ):=/jJ(6i-<Oi))\ (22) 

\* =1 / Weyl 

a complete set of variables results. Here, the subscript "Weyl" in the definition of the 
moments indicates that we are ordering all operator products totally symmetrically. Any 
of the basic operators, f3+, p+, a and p a , can appear as the Oi in the basic moments. (To 
match with standard notation, we will write A(0 2 ) = (AO) 2 for fluctuations.) 

Expectation values of basic operators together with the moments can be used to char- 
acterize an arbitrary state (pure or mixed); they can be used as coordinates on the state 
space. Moreover, the commutator of basic operators endows the state space with a Poisson 
structure, defined by 

{(i),(5)}:=M^!l (23) 

for any operators A and B. By linearity and the Leibniz rule, this defines Poisson brackets 
between all the moments and expectation values. In this way, the quantum phase space 
is defined. It is not a linear space since there are restrictions for the moments, most 
importantly the uncertainty relations such as 

(A«) 2 (Ap Q ) 2 -(A(ap Q ) 2 >^. (24) 

On this kinematical quantum phase space, the constraint C must be imposed. For a 
constraint operator polynomial in the basic operators, 

C pol :=((p^i-(pd))(7) (25) 
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defines ajhnction on the state space, expandable in the moments, for any arbitrary poly- 
nomial pol of the basic operators. This infinite set of functional satisfies two important 
properties: (i) all these functions vanish on physical states annihilated by C, and (ii) they 
form a first class set in the sense of classical constraint analysis, i.e. {C po \, C po y} vanishes 
on the subset of the phase space where all C po \ vanish. The quantum constraint C can 
thus be implemented on the space of moments by imposing the infinite set {C po i} as con- 
straints as one would do it on a classical phase space. We have to find the submanifold on 
which all constraints vanish, and factor out the flow generated by them. If this reduction 
is completed, we obtain the physical quantum phase space and can look for solutions of 
observables. 

This procedure has several advantages [6]. As already mentioned, it is completely 
representation independent and instead focuses on algebraic aspects of a quantum system. 
As a consequence, implementing constraint operators with zero in their continuous spectra 
is no different from implementing those with zero in their discrete spectrum. Any difficulties 
in finding physical inner products can be avoided, for the physical normalization arises 
automatically when the constraints are solved for moments. 

Once we try to find specific solutions, there are of course practical difficulties. We are 
dealing with infinitely many constraints on an infinite-dimensional phase space. Sometimes, 
this set of equations decouples to finitely many ones, but this happens only in rare solvable 
cases. In more general systems, we must use approximations to reduce the set to a finite 
one of relevant equations, with a semiclassical approximation as the main example. Here, 
we look for solutions whose moments satisfy a certain hierarchy, higher moments in the 
semiclassical case being suppressed by powers of h, A(Oi...0„) oc h n / 2 . To any given 
order in h, only a finite number of moments need be considered, subject to a finite number 
of non-trivial constraints. 

In our model, to second order in moments we have the effective constraints 



c 


= (p + ) 2 -(p a ) 2 - 


- (a) 2 + (Aj9 + ) 2 - {Ap a f - {Aaf = 


(26) 




= 2(p + )A((3 +P+ ] 


) + ih{p + ) - 2(p a )A((3 +Pa ) - 2(a)A((3 + a) = 


(27) 


c p+ 


= 2(p + )(A P+ ) 2 - 


- 2(p a )A( P+Pa ) - 2(a)A(p + a) = 


(28) 




= 2(p + )A( P+ a) 


- 2(p a )A(ap a ) - ih(p a ) - 2(a)(Aa) 2 = 


(29) 




= 2{p + )A{p + p a ) 


- 2(p a )(A Pa ) 2 - 2(a)A(a Pa ) + ih(a) = . 


(30) 



The reduction has been performed in [7], with the result that the physical state space is 
equivalent to a deparametrized quantum system with constraint Cq = (p + ) ±Hq with the 
reduced Hamiltonian 

n rp v—r^ fi , («) 2 (Ap a ) 2 - 2{a){p a )A{ap a ) + {p a ) 2 {Aa) 2 \ 

= V (Pa> 2 + + 2{{pa)2 + (&)2)2 J (31) 

+0 {(Ap a ) 4 ) + O ((Aa) 4 ) + O {(A(ap a )) 2 ) . (32) 

A similar setup can be used in BRST quantizations, where the BRST charge would replace the con- 
straint operators [TU] , 
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To leading order, we reproduce the classical Hamiltonian, but there are corrections from 
quantum fluctuations and correlations coupling to the expectation values. Moreover, from 
the Poisson brackets between moments we obtain Hamiltonian equations of motion telling 
us how a state spreads or is being squeezed. We have 

<*(<*) (Pa) / 33 x 

y/{p a ) 2 + <«> 2 

, (p a )(Aa) 2 (2(a) 2 - (p a ) 2 ) + (a)A(ap a ) (4(p a ) 2 - 2(a) 2 ) - 3(p a ) (u) 2 (Ap a ) 2 



2«p a > 2 + (a> 2 )^ 

d(p a ) _ -(«) 



(34) 



, 3(«)(p Q ) 2 (Aa) 2 - (p Q )A(«j9 Q ) (4(a) 2 - 2<p a ) 2 ) - (d)(Aj9 Q ) 2 (2(p a ) 2 - (a) 2 ) 



d(Aa) 2 2(a) 2 A(ap a ) - 2(a)(p a )(Aa) 2 

((p a ) 2 + («) 2 ) f 
d(Ap Q ) 2 2(&)(p a )(Ap a ) 2 -2(a) 2 A(ap a ) 



dA(a Pa ) = (a) 2 (Ap a ) 2 - (p a ) 2 (Aa) 2 
((j5 a ) 2 + («) 2 ) f 



(35) 
(36) 
(37) 



Solutions (a)(/3 + ), (p a )(/3+), (Aa) 2 (/3 + ), (Ap a ) 2 (/? + ) and A(ap a )(/3 + ) provide physical 
observables, corresponding to expectation values and moments in physical states, telling 
us how a state moves and spreads. 



3 Comparison 

Initially, the kinematical coherent state with expansion coefficients given by (|2"T|) yields the 
expectation values (d)| t= o = Re(z) and (p a )\t=o = Im(z). The second order moments then 
saturate the uncertainty relation and their values are 

(Aa) 2 | t=0 = ~ (Ap a ) 2 \ t=0 = ^ A(ap a )| t=0 = 0. 

For a concrete comparison we select a state that is initially peaked about (a) = 
and (p a ) = 0, so that z = a . In order for the state to be semiclassical, we need a to be 
"significantly larger" than y/h. We make a concrete choice ao = IOTP. In Fig. [TJ we plot 
alongside each other the classical and two quantum-corrected trajectories of the system 
starting from the above initial values. The two corrected trajectories were calculated using 
two different methods: the kinematical coherent state of Section [2721 and effective equations 
truncated at order h, Eqs. (l3~3l-(!37l) of Section 12.31 
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The two quantum trajectories agree very well for much of the evolution shown. In 
Figs. [2HH we plot the quantum evolution of the second order moments generated by a and 
p a . From Fig. [3] in particular it is clear that the semiclassical approximation breaks down 
somewhere between /3 + = 2a$ and j3 + = 3«o as Ap a is no longer "much smaller" than 
qjo- Up until that point both methods for quantum evolution are in close agreement in 
describing not only the trajectory in the a — p a space but also the evolution of the second 
order moments themselves. Since semiclassicality was used to obtain the truncated system 
of equations, there is no reason to expect it to be accurate beyond f3 + = 3ao. 



Aa 




, , , , 1 , , , , 1 , , , , 1 , , , , 1 — yS + 

a o 2a 3a 4a 

Figure 2: Coherent state (solid) and effective (dashed) evolution of the second order mo- 
ment Aa = y/ (Aa) 2 in units of «o- 



The two methods agree excellently where the domains of their applicability overlap. 
Both methods have their own strengths and weaknesses. The weakness of the semiclassi- 
cally truncated effective equations should by now stand out — in some systems semiclas- 
sicality eventually breaks down and moments of high order dominate. Direct evaluation 
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Figure 4: Coherent state (solid) and effective (dashed) evolution of the second order mo- 
ment A(apa) in units of ckq. 
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on the states does not rely on this approximation and remains a valid method for com- 
puting the wave function at all times. The shortcomings of the latter technique are less 
obvious and are of practical nature. In order to apply the method one needs, first of all, 
to decompose the initial state as a sum of the eigenstates of the Hamiltonian, which for 
an arbitrary state and a typical Hamiltonian can be complicated. During evolution, each 
of the eigenstates acquires a phase factor and they need to be re-summed to compute the 
wave function at a later time. For an arbitrary state, the sum may converge very slowly 
requiring one to sum over a very large number of eigenstates to obtain an accurate de- 
scription of the wave function. Finally, for expectation values and moments of observables 
further integrations must be done. In systems of several degrees of freedom, this will add 
considerably to computation times. 

If one is to make robust predictions, a range of initially semiclassical states with a variety 
of initial values of moments should be considered — the procedure is very complicated to 
implement using state decomposition but amounts to nothing more than simply changing 
the initial conditions in the case of the effective equations. In the subsections that follow 
we use the methods separately and exploit their individual strengths. 

3.1 Long-term behavior of the state 

Knowing the state exactly, allows us to plot the magnitude of the wave function and make 
precise long-term predictions. From Fig. Owe see that after j3 + ~ 10«o the state becomes 
highly quantum, spread out over an entire orbit. Expectation values can no longer be 
interpreted as the most probable outcome of a measurement. 



In Figs. EHH] we look at the long term behavior of the leading order quantum degrees 
of freedom. The magnitudes of (Aa) 2 and (Ap a ) 2 rise rapidly until around [3+ = 40ao. 
The state remains stable during 40ao < f3 + < 200«o after which all second order moments 
begin to oscillate with an increasing amplitude. The amplitude reaches its peak around 
P + = 300ao; thereafter the moments keep oscillating with a stable amplitude for as long as 
the evolution has been traced, up to around /3 + = 10 5 ao- As can be seen from the example 
of ((& — (a)) 3 ) in Fig. O third order moments follow a similar pattern. 
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Figure 5: Square-magnitude of the coherent state wave function evolved for < (3 + < 15a . 
Aa 




50 100 150 200 250 300 



Figure 6: Coherent state evolution of the second order moment Aa = a/ (Aa) 2 in units of 
a . 
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Figure 7: Coherent state (solid) and effective (dashed) evolution of the second order mo- 
ment Ap a = a/ (Ap a ) 2 in units of a . 
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Figure 8: Coherent state (solid) and effective (dashed) evolution of the second order mo- 
ment A(ap a ) in units of a . 
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Figure 9: Coherent state evolution of the third order moment ((& — (a)) 3 ) in units of a®. 

As expected, semiclassicality eventually breaks down and quantum fluctuations become 
large; this trend is illustrated by the evolving uncertainty measure (Aa) 2 (Ap a ) 2 — (A(ap a )) 2 
(bounded below by the uncertainty relation) shown in Fig. [TD1 After the initial increase, 
there is a period of approximate stability; then the leading order fluctuations start to 
oscillate with large amplitudes. Even though some moments return to small values during 
these oscillations, semiclassicality is not regained as shown be the long-term behavior of 
the uncertainties in Fig. [TTJ Details of this behavior, found numerically, appears rather 
characteristic, but it is difficult to find an explanation based on the underlying equations. 



3.2 Short-term evolution with varying initial conditions 

We now evolve effectively starting from the same initial expectation values, but varying the 
initial values of second order moments. By the specific choices, some sets of moments used 
here no longer saturate the uncertainty relations, and some have non- vanishing correlations. 
They cannot correspond to Gaussian states, and so their initial configurations would be 
much more difficult to implement using wave functions of physical states. The results are 
plotted in Fig. [TJJ In Figs. fTSHTol we plot the corresponding effective evolution of moments 
for each set of the initial conditions, where the thin horizontal line approximately indicates 
the threshold of the semiclassical approximation for the second order moments at .04«q. 
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Figure 10: Short-term coherent state evolution of the "quantum uncertainty" defined as 
(Aa) 2 (A Pa ) 2 - {A{ap a )f in units of f . 
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Figure 11: Long-term coherent state evolution of the "quantum uncertainty" defined as 
(Aa) 2 (Ap a ) 2 - (A(ap a )) 2 in units of f . 
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Figure 12: Effective evolution of semiclassical states with different initial values of moments 
for < f3 + < 3.5ao. The initial values of moments are as follows. Dashed line — the original 
coherent state: (Aa) 2 = (Ap a ) 2 = .005a 2 ,, A(ap a ) = 0. Dotted line: (Aa) 2 = .025a 2 ,, 
(Ap a ) 2 = .001a 2 ,, A(ap a ) = 0. Thin solid line: (Aa) 2 = .001a 2 ,, (Ap a ) 2 = .025a 2 ,, 
A(ap a ) = 0. Thick solid line:(Aa) 2 = .008a 2 ,, (Ap a ) 2 = .008a 2 , A(ap a ) = .005a 2 ,. 
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Larger (Aa) 2 results in a larger deviation from the classical behavior and a faster 
breakdown of the semiclassical approximation. The above effect is much less sensitive to 
the momentum spread (Ap a ) 2 . This disparity between the effects of the two spreads may 
seem surprising given the symmetry between a and p a in the expression for the physical 
Hamiltonian. This symmetry, however is broken by the initial state we have chosen, which 
is peaked about (p a ) = and (a) = ao, so that the spread in a, produces a larger spread 
in the energy than the spread in p a . 
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Figure 13: Effective evolution of second order moments (Aa) 2 (dashed), (Ap a ) 2 (solid) 
and A(ap a ) (dotted), which initially take the values .005a 2 ,, .005a 2 , and respectively. 



4 Effective evolution of a radiation-filled universe 

In this section we use the effective equations on their own to analyze quantum corrections 
to the dynamics of a radiation-filled Bianchi I universe briefly mentioned in Section [3 This 
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Figure 14: Effective evolution of second order moments (Aa) 2 (dashed), (Ap a ) 2 (solid) 
and A(ap a ) (dotted), which initially take the values .025a: 2 ,, .OOla: 2 , and respectively. 
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Figure 15: Effective evolution of second order moments (Act) 2 (dashed), (Ap a ) 2 (solid) 
and A(ap a ) (dotted), which initially take the values .OOla: 2 ,, .025a 2 , and respectively. 
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Figure 16: Effective evolution of second order moments (Aa) 2 (dashed), (Ap a ) 2 (solid) and 
A(ap a ) (dotted), which initially take the values .008a 2 ,, .008a 2 , and .005a 2 , respectively. 

model has a more realistic matter content than the one analyzed above. We recall that the 
energy density of radiation has the form p oc a -4 = e~ Aa , which results in the Hamiltonian 
constraint C = p 2 + — p 2 a + e 2c V. The constraint condition may be implemented effectively 
in fashion a very similar to the way it was done in Section 12.31 however, the physical inner 
product treatment would require a detailed knowledge of the (now continuous) spectrum of 
an operator that is completely different from H that was used in (fT5j) and is not straight- 
forward to obtain. Even if we could determine energy eigenstates, expanding Gaussians or 
other general semiclassical states in this basis would be challenging. For this reason we do 
not attempt to implement the constraint of the radiation-filled model on the kinematical 
Hilbert space and restrict our analysis to the effective procedure, demonstrating its wide 
applicability. 

4.1 Classical behavior 

We deparametrize the constraint exactly as in Section I2.1[ selecting (3 + as time. The 
dynamics on a and p a is then generated by the Hamiltonian H = \/p 2 a — e 2c V and results 
in the equations of motion da/d(3 + = p a / H and dp a /d(3 + = e 2c V/iJ. Noting that H is once 
again a constant of motion we can immediately infer the classical phase-space trajectories, 
which are of the form p a = ±V 'e 2oL v + const. They split the a — p a space into three regions 
as illustrated in Fig. [T71 There are no classical orbits in Region 3, as H becomes complex. 
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Figure 17: Disjoint regions of classical solutions. Region 1: trajectories have the shape 
p a = V 'e 2a v + const, they correspond to expanding universes as a increases with /3 + . 
Region 2: trajectories have the shape p a = —\/ e 2a u + const, they correspond to contracting 
universes as a decreases with f3+. Region 3: no classical solutions. 
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The explicit solution to the equations of motion for evolution in terms of (3 + is given by 

a{(3 + ) =/3 + -log(2(l + Ae 2/3 +)) +B (38) 

where the integration constants A and B can be related to the initial values «o := a(0) 
and p a0 : = p a (0) via 

B = log ( 4 + IS {J + «o 

In terms of A and S, the constant of motion is 

which requires A < in order for H to be real. Using ( 1381) . p a (/3 + ) may be recovered from 
the equations of motion (/?+) = Hda/df3 + . 

For the orbits in Region 1, a and p a reach infinity at a finite positive value of the 
evolution parameter, namely at (3 + = «o — log(p a o — H). Explicit integration of the 
expression (flTI) for the proper time using the above solution for a(/3 + ) shows that this 
value of (3 + is reached in infinite proper time. In other words, the expansion, as expected 
for the radiation filled universe, takes infinite proper time, but anisotropy asymptotically 
approaches a maximum value. For the orbits in Region 2 one can obtain a similar result by 
tracing evolution backwards in time; in this case a reaches infinity and p a reaches negative 
infinity at (3 + = cto — log(|p a o| + H), which is always negative in that region. 

Following those orbits forward in time, one finds that the collapse happens only as (3 + 
goes to infinity. Once again one can use ( TTTT) to convert this to a proper time interval, with 
the result that, as one would expect, the collapse takes a finite amount of proper time, 

n 3 p 2B /I IS 

At := T (P + = oo) - r(J3 + = 0) - ^ 1 '' !| 




32^% \ {\A\ -l) 2 
In this model, with positive energy, the singularity is certainly not resolved. 

4.2 Effective constraints 

Following the effective procedure for solving constraints outlined in Section I2T31 we find the 
constraint functions truncated at second order: 

C = (p + ) 2 -(p a ) 2 + e 2{&) v + (Ap + ) 2 -(A Pa ) 2 + 2e 2{&) v(Aa) 2 = (39) 

C p+ = 2(p + )A(P + p + ) +i h(p + )-2(p a )A(f3 + p a ) + 2e 2 ^vA([3 + a)=0 (40) 

C p+ = 2(p + )(Ap + ) 2 -2(p a )A(p +Pa ) + 2e 2 ^uA(p + a) = (41) 

C a = 2(p + )A(p + a)-2(p a )A(ap a )-th(p a ) + 2e 2 ^u(Aa) 2 = (42) 

C Pa = 2(p + )A(p + p a )-2(p a )(Ap a ) 2 + 2e 2{&) vA(ap a ) -thesis = 0. (43) 
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These can be solved and gauge fixed following the same steps as for the model of Section [231 
with the result that evolution in (3 + on the expectation values and moments of d and p a 
is generated by the quantum Hamiltonian 

_ J, t y. W x , c2 (,y -( W + 2(p a )A(a Pa ) + (e^u - 2(p a ) 2 )(Aa) 2 \ 

+0 ((A Pa ) 4 ) + O ((Aa) 4 ) + O ((A(a Pa )) 2 ) . (44) 

The equations of motion are obtained by taking the quantum Poisson bracket between 
quantum variables and the quantum Hamiltonian. 

d(d> (Pa) (45) 



, C 2(*)„ (P°> ((&»>' + ( A «) 2 ~ A (^") ( 2 ^) 2 + ~ i(P«)(APa) S 



d/3+ V^) 2 " e2<A> ^ 

: )2 + l e 2(d>^ (Ao 

d^> = -(d) (46) 

d/3+ v /(p a )2 _ e 2< (i > z/ 

, c 2(^, (l e4< ^ 2 - (Pa) 2 e 2 < & ^ + 2(p Q ) 4 )(AoQ 2 - 2A(a Pa ) (\e 2 ^v + (p Q ) 2 ) (p a ) 

{(Pa) 2 -e 2 (°)z/)5 

e ^;ApJ 2 (6U 2 + ±e 2 < a M 



;2(A) JA Po ) 2 ((p a ) 2 + le 2 ^) 

2((p Q ) 2 -e 2 ^)z/) f 
d(Aa) 2 2e 2 ^z/((p a )(Aa) 2 - A(qp Q )) 

dp+ ((p a ) 2 -e 2 ( & )is)^ 
d(Ap a ) 2 2e 2 ^v ((2(p a ) 2 - e 2 ^v) A(ap a ) - (p a )(Ap a ) 2 ) 
d(5 + ^ a )2_ e 2{&) v y 

dA{ap a ) e 2 ^u ((2(p a ) 2 - e 2 ^v) {Aaf - (Ap a ) 2 ) 



df3+ «p a ) 2 -e 2 <*>z/)* 

As before, these equations reduce to the classical equations of motion at zeroth order in 
moments and are straightforward to evolve numerically. 

4.3 Numerical evolution 

In this section we take a classical phase space trajectory from each of the Regions 1 and 2 
and compare the effective trajectories for a variety of semiclassical states initially peaked 
about these classical solutions. Effective and classical trajectories of the expanding universe 
are plotted in Fig. [181 Clearly, the significance of quantum back-reaction can be seen well 
before the approximation breaks down. 



(47) 
(48) 
(49) 
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Figure 18: Phase space trajectories of an expanding universe classical (dotted) and effective 
with different initial values for second order moments: dashed line — (Aa) 2 u = .OOlHj 2 , 
(Ap a ) 2 = M5H 2 , A(ap a )^/V = 0; thin solid line— (Aafu = .025# 2 , (Ap a ) 2 = .001H 2 , 
A(ap a )^ = 0; thick solid line— {Aa) 2 v = .008H 2 , (Ap a ) 2 = .008H 2 , A(ap a )^ = 
.005Hq. Solutions were evolved for < (3 + < .2; the vertical line indicates the breakdown 
of the semiclassical approximation based on the size of second order moments. 
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The corresponding evolution of the leading order moments starting from different initial 
values is plotted in Figs. [1911211 where the horizontal line, as before, indicates an approxi- 
mate threshold of the semiclassical approximation for the second order moments. 



Moments 




Figure 19: Evolution of second order moments in an expanding universe. (Aa) 2 i/ (dashed), 
{Ap a f (solid), A(ap a )^/v (dotted), with initial values: (Aa) 2 v = .008H 2 , {Ap a f = 
.008# 2 , A(ap a )^- 



.005H 2 . 



Effective and classical trajectories of the contracting universe are plotted in Fig. [221 The 
corresponding evolution of the leading order moments starting from different initial values 
is plotted in Figs. [2311251 where the horizontal line, once again, indicates an approximate 
threshold of the semiclassical approximation for the second order moments. 
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Figure 21: Evolution of second order moments in an expanding universe. (Aa) 2 u (dashed), 
(Ap a ) 2 (solid), A(ap a )y/u (dotted), with initial values: (Aa) 2 v = .001# 2 , (Ap a ) 2 = 
.025H 2 , A(a Pa )^=0. 
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Figure 22: Phase space trajectories of a contracting universe classical (thin solid line) 
and effective with different initial values for second order moments: dashed line— 
(Aa) 2 u = .001F 2 , (Ap a ) 2 = .025# 2 , A(ap a )^ = 0; dotted line— (Aafv = .025# 2 , 
(Ap Q ) 2 = .001F 2 , A(ap a )y/v = 0; thick solid line— {Aafv = .008F 2 , (Ap a ) 2 = .008H 2 , 
A(ap Q ,)y / z/ = .005Hq. Solutions were evolved for < f3 + < .07; the vertical line indicates 
the breakdown of the semiclassical approximation (coming from the right along the a- axis). 
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Figure 25: Evolution of second order moments in a contracting universe. (Aa) 2 z/ (dashed), 
(Ap a ) 2 (solid), A(ap a )^/v (dotted), with initial values: (Aa) 2 v = .001if£, (Ap a ) 2 = 
.025H 2 , A(ap a )y/V=0. 

5 Conclusions 

Most quantum systems of physical relevance can only be analyzed by perturbation meth- 
ods. Quantum gravity and quantum cosmology cannot be considered exceptions. As 
demonstrated by the examples of this article, canonical effective equations, based on the 
back-reaction of moments of a state on its expectation values, are of wide applicability, 
capture quantum effects reliably, and approximate the full quantum dynamics in a self- 
consistent way. They are tractable even in situations where semiclassical wave functions in 
physical Hilbert spaces would be too complicated to be constructed — cases which abound 
in quantum cosmology and quantum gravity. 
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